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Abstract 

Action potential duration (APD) restitution, which relates APD to the preceding diastolic in- 
terval (DI), is a useful tool for predicting the onset of abnormal cardiac rhythms. However, it is 
known that different pacing protocols lead to different APD restitution curves (RCs). This phe- 
nomenon, known as APD rate-dependence, is a consequence of memory in the tissue. In addition 
to APD restitution, conduction velocity restitution also plays an important role in the spatiotem- 
poral dynamics of cardiac tissue. We present new results concerning rate-dependent restitution in 
the velocity of propagating action potentials in a one-dimensional fiber. Our numerical simulations 
show that, independent of the amount of memory in the tissue, waveback velocity exhibits pro- 
nounced rate-dependence and the wavefront velocity does not. Moreover, the discrepancy between 
waveback velocity RCs is most significant for small DI. We provide an analytical explanation of 
these results, using a system of coupled maps to relate the wavefront and waveback velocities. Our 
calculations show that waveback velocity rate-dependence is due to APD restitution, not memory. 
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I. INTRODUCTION 



When a cardiac cell is depolarized by an electrical stimulus, it exhibits a prolonged 
elevation of transmembrane potential known as an action potential. We define the action 
potential duration (APD) as the time required for the cell to achieve 80% repolarization 
following a depolarizing stimulus. The refractory period between the end of an action 
potential the application of a subsequent stimulus is called the diastolic interval (DI). It is 
known that APD restitution (the dependence of the APD on preceding DI) is of fundamental 

tin 

importance in paced cardiac dynamics. In particular, studies |2J show that the slope of the 
APD restitution curve (RC) is linked to the onset of alternans, an abnormal cardiac rhythm 
characterized by long-short variation of APD, which may lead to ventricular fibrillation and 
sudden cardiac death y, ]4 , M . 

Experts* flflfi and ana. ytl eal Q inventions have shown that different pacing 



protocols lead to different APD RCs, a phenomenon known as APD rate-dependence. Several 
studies indicate that the origin of APD rate-dependence is the presence of memory 

in cardiac tissue. That is, APD depends not only upon the preceding DI but also on the 
previous history of paced cardiac tissue. Memory appears to be a generic feature of cardiac 
muscle since it has been reported in humans [^] and various animals jl, 13, 14, 



Testing for rate-dependence involves the use of multiple pacing protocols and comparison 
of the resulting RCs. Two of the most commonly used pacing schemes are the dynamic and 
S1-S2 pacing protocols. Under the dynamic (steady-state) protocol, pacing is performed at 
a constant basic cycle length B until steady-state is reached (no beat-to-beat variation in 
APD or DI). After recording the steady-state DI-APD pair, B is changed by an amount A 
and the process is repeated. The dynamic RC is constructed by plotting all steady-state 
DI-APD pairs obtained from the dynamic pacing protocol over a range of B values. The Sl- 
S2 (standard) protocol also begins with pacing at a fixed basic cycle length B (SI interval) 
until steady-state is reached. Then, an S2 stimulus is applied at an interval B\ after the 
final SI stimulus. Setting 5 = B\ — B , the S1-S2 RC is obtained by plotting the APD 
following the S2 stimulus versus the preceding DI for different values of 5. Note that there 
is only one dynamic RC, whereas each different SI pacing interval can yield a distinct S1-S2 
RC. For the purposes of this paper, we obtain only local S1-S2 RCs (\5\ small relative to 
SI) for different values of the SI interval. In particular, following ^| we apply one short 
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FIG. 1: Typical APD restitution curves obtained using dynamic and S1-S2 pacing protocols (A = 
40 ms, 5 = ±20 ms) for two different ionic membrane models, (a) A two-current model, (b) A 
three-current model. The dynamic RC (solid) and local S1-S2 RCs (dashed) are shown for different 
values of SI. 

(Bi = B — 5) and one long (B\ = B + 5) perturbation at each different value of B. 

The connection between rate-dependent APD restitution and memory is illustrated in 
Fig. H which shows APD RCs obtained from numerical simulations using two different ionic 
models (see Section II) of the cell membrane. Figure ^ is generated using a two-current 



model 



» (see Se 

HQ 



with no memory: the dynamic and S1-S2 RCs are indistinguishable. Figure 
is generated using a three-current ionic model 3, Q with some memory. One can see 
from Fig. [I]p that segments of S1-S2 RCs (dashed curves) do not coincide with the dynamic 
RC (solid curve) for small DI values. The splitting between the dynamic and S1-S2 RCs is 
the manifestation of APD rate-dependence. 

In the case of a network of electrically-coupled cardiac cells, conduction-velocity resti 
tution also plays an important role in dynamics of the spatially extended system 
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Indeed, when a sequence of propagating pulses is produced, both APD and the propagation 



speed of a pulse are influenced by the preceding pulse. Conduction-velocity restitution is 
analogous to APD restitution in that it relates the speed of an action potential at a given 



site to the DI at that site. Several authors 
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24j have noted that abrupt changes in 



the pacing rate lead to discrepancies between the wavefront and waveback velocities. For 
this reason, we will always distinguish between wavefront (or activation front) velocity and 
waveback (or recovery front) velocity of propagating action potentials. 

In this paper, we investigate rate-dependence of wavefront and waveback velocities of 
propagating action potentials in a one-dimensional fiber of cardiac cells. Using numerical 
simulations of different ionic membrane models (with and without memory) we demonstrate 
that the waveback velocity exhibits pronounced rate-dependence and the wavefront velocity 
does not (Section II). We derive an analytical relationship between wavefront and waveback 
velocities by modeling cardiac dynamics using a system of coupled maps with an arbitrary 
amount of memory. We show that APD restitution, not memory, leads to waveback velocity 
rate-dependence (Section III). We provide conclusions and discussion in Section IV. An 
Appendix on the two-current ionic model is included for reference. 



II. RATE-DEPENDENT VELOCITY: NUMERICAL RESULTS 

Typically, the cardiac action potential is modeled by considering ionic currents that flow 
across the cell membrane via ion channels. The rate-of-change of the transmembrane voltage 
is obtained by summing all ionic currents and dividing by the membrane capacitance. The 
ion channels act as gates that regulate the permeabilities of ions, most notably sodium, 
potassium and calcium. Hence, ionic models are presented as systems of ordinary differential 
equations that govern transmembrane voltage and gate variables. 

In the case of a one-dimensional fiber, electrical coupling can be modeled by the inclusion 
of a diffusion term. The result is a reaction-diffusion partial differential equation known as 
the cable equation: 

dv _ d 2 v J total 0<T</ m 

where v denotes transmembrane voltage, x measures distance from the stimulus site, C m is 
membrane capacitance, k is a diffusion coefficient, and /total is the sum of all ionic currents. 
The number of currents varies depending upon the complexity of the ionic model. The 
diffusion coefficient incorporates membrane capacitance, cell surface-to-volume ratio, and 
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longitudinal resistivity of cardiac muscle tissue. In all of our numerical simulations C m = 
1 /xF cm -2 and k = 0.001 cm 2 ms" 1 . Neumann boundary conditions v x (0,t) = v x (L,t) = 
are imposed at both ends of the cable. 

To investigate rate-dependent propagation, we perform numerical simulations of Eq. (JTJ). 
We apply both dynamic and S1-S2 pacing protocols at one end of a cable, and measure the 
wavefront and waveback velocities of each propagating pulse. By analogy with APD rate- 
dependence, velocity rate-dependence means that different pacing protocols lead to different 
velocity RCs. 

Since memory is responsible for APD rate-dependence, it is natural to hypothesize that 
memory also leads to wavefront and waveback velocity rate-dependence. Consequently, we 
use two different ionic models in our numerical simulations: a two-current ionic model ^| 
with no memory and a three-current ionic model [sij with some memory. 

The details of the numerical experiments are as follows. Using a cable of length L = 10cm, 
we solve Eq. Q numerically with an operator-splitting method. Stimuli are applied over 
a 1 mm region at the proximal (x = 0) end of the fiber using both the dynamic and Sl- 
S2 protocols described in the Introduction. In all simulations, we use A = 40 ms and 
5 = ±20 ms. Pacing results in a train of pulses that propagate left-to-right in the fiber. 
Measurements of DI, APD, wavefront speed, and waveback speed are taken at x = 2.5 cm. 
The position of a pulse wavefront is defined as the x value for which the transmembrane 
voltage is — 60mV and dv/dx < 0. Likewise, waveback position is defined as the x value 
at which v = — 60mV and dv/dx > 0. Linear interpolation is used to improve tracking of 
wavefront and waveback positions. Speeds are then computed by recording the time required 
for wavefronts and wavebacks to traverse a 1-mm-wide interval centered at x = 2.5 cm. For 
illustration purposes, Fig. |21 shows a projection of a steady-state solution of Eq. (with 
two-current ionic model) onto the xt plane. Different shades of grey correspond to different 
transmembrane voltages, with black corresponding to the rest potential. Note that, in 
steady-state, projecting the wavefronts and wavebacks onto the xt plane forms a sequence 
of parallel lines. 

Results of numerical simulations of Eq. Q with the two-current ionic model (no memory) 
are presented in Fig. El which shows wavefront (Fig. |3^) and waveback (Fig. Eb) velocity 
RCs. One can see from FigEK that the wavefront velocity RCs resulting from different pacing 
protocols are indistinguishable. Thus, there is no significant rate-dependence if velocities 
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FIG. 2: Steady-state space-time response obtained by numerical simulation of Eq. (^Q) with the 
two-current ionic model. 

are measured at the wavefront. However, one can see from FigEb that segments of S1-S2 
waveback velocity RCs (dashed curves) do not coincide with the dynamic waveback velocity 
RC (solid curve). As in the case of APD rate-dependence (see Fig. [TJ), the splitting between 
waveback velocity RCs is more pronounced for small values of DI. 

Wavefront and waveback velocity RCs obtained from numerical simulations of a cable 
with the three-current ionic model (that has some memory) are presented in Fig. The 
results are qualitatively similar to the two-current model results shown in Fig. El 

There are two important points that we wish to emphasize. First, rate-dependent wave- 
back velocity restitution does not depend upon the presence of memory in the tissue, as 
evidenced by our two-current model simulations. Second, rate-dependent waveback veloc- 
ity is more pronounced for small values of DI. In what follows, we provide an analytical 
explanation of these findings. 



III. RATE-DEPENDENT VELOCITY: ANALYTICAL RESULTS 



Instead of considering the systems of ODEs that constitute ionic membrane models, many 
authors employ mapping models that describe APD as a function of the previous DI and 

25^ . To our knowledge, Nolasco and Dahlen were the first to 
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propose a simple mapping model of the form 



A 



n+l 



F{D n ) 



(2) 



to describe cardiac dynamics. Here, A n and D n denote the n th APD and DI values, re- 
spectively. In general, the number of arguments of F determines how much memory is 
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FIG. 3: Wavefront and waveback velocity RCs obtained from numerical simulation of Eq. Q with 
two-current ionic model (A = 40 ms, 5 = ±20 ms). Velocities were measured at x = 2.5 cm. (a) 
Wavefront velocity RCs. (b) Waveback velocity RCs. The dynamic RC is solid and the local S1-S2 
RCs are dashed. 



present. 

For certain ionic models, such as the two and three-current models, it is possible to derive 
mappings by analyzing the ODEs. As demonstrated in a mapping of the form (J2J) can 
be derived directly from the two-current model equations (see Appendix). The mapping 
model (j2J) has no memory, and all APD RCs coincide as in Fig. QJi. 

It was shown in j3 that the three-current model leads to a mapping with two arguments: 



A n+1 = F(A n ,D n ). 



(3) 



This mapping model has some memory, and the dynamic and S1-S2 RCs are different as in 
Fig. do. 

In order to explain differences between wavefront and waveback velocity RCs for different 
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FIG. 4: Wavefront and waveback velocity RCs obtained from numerical simulation of Eq. Q with 
three-current ionic model (A = 40 ms, S = ±20 ms). Velocities were measured at x = 2.5 cm. (a) 
Wavefront velocity RCs. (b) Waveback velocity RCs. The dynamic RC is solid and the local S1-S2 
RCs are dashed. 

pacing protocols, we approximate the dynamics of Eq. (JTJ) with a system of coupled maps. 
We follow the approach described in 26[, which allows us to derive a relationship between 
wavefront and waveback velocities. We analyze the dynamic and S1-S2 pacing protocols 
separately since the pacing protocol determines the boundary conditions for Eq. (|T)). 



A. Dynamic pacing protocol 

Under the dynamic pacing protocol, pacing is performed at a constant basic cycle length, 
B, at x = until steady-state is reached. In what follows, we assume that a 1:1 steady-state 
response results from dynamic pacing. That is, every stimulus produces an action potential 
and there is no beat-to-beat variation in APD or DL A schematic representation of steady- 
state behavior is shown in Fig. El which shows the projection of a particular level set of the 
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FIG. 5: Schematic representation of a 1:1 steady-state. Here, 4> n denotes the n activation front 
and (5 n is the corresponding waveback. Solid lines represent wavefronts and dashed lines represent 
wavebacks. 

surface in Fig. El onto the xt plane. The lines in Fig. are identified with the sequence of 
wavefronts and wavebacks. We define <f> n (x) ((3 n (x)) as the time at which the n th wavefront 
(waveback) reaches x. Note that 4> n {x) and f3 n (x) are parallel lines in the xt plane if a 1:1 
steady-state is reached. The cycle length is defined as 

CL n (x) = A n (x) + D n (x) = 4> n +i(x) - 4> n (x). (4) 

We remark that CL n (0) = B for all n, and CL n (x) = B for all x if steady-state is reached. 

Let us assume that, at each x along the fiber, APD can be represented as a function of 
an arbitrary number of preceding APDs and DIs in a form 

A n+l {x) = F{A n {x),B n {x)), (5) 

i.e. an arbitrary amount of memory is included. Here, 

A n (a;) = (A n (x),A n -i(x), ...,A n+1 - m (x)), (6) 
D„(x) = (D n (x),D n - 1 (x),...,D n+1 - k (x)), 

and m > and k > 1 are integers characterizing how many preceding states are taken into 
account in the mapping model. Since many previous states are involved, Eq. ((HJ) makes 
sense only for m,k < n. Note that m — 0, k — 1 corresponds to the simplest mapping 
model Eq. (J2J), the case of no memory. The case m — 1, k = 1 corresponds to a mapping of 
the form of Eq. Q with some memory. 

9 



Let us also assume that the velocity of the (n + l) st wavefront, Cf ront (A n (x), D n (x)), 
depends upon preceding (local) APD and DI values. This velocity is computed by inverting 
the slope of 4> n (x): 

= 7T n n r ft = G ^ x ^ W 

dx Cfr OIlt (A n (a;),D„(a;)) 

When steady-state is reached, the vectors A n (x) and D n (x) are constant: 

A n (x) = A* = (A*,A*,...A*), (8) 
B n (x) = D* = (D*, D*, . . . D*). 

Thus, plotting Cf ront (A*, D*) versus D* , we obtain a point on the dynamic wavefront velocity 
RC. The curves f3 n+ x(x) and <f) n+ x(x) have the same slope since they are parallel at steady- 
state: (3 n+1 (x) = <f> n+ i(x) + A*. Therefore, the dynamic waveback and wavefront velocity 
RCs are identical. Hence, from now on we refer to the dynamic velocity RC and use the 
notation Cd yn = Cd yn (-D*)- Since propagation speeds typically increase when more recovery 
is allowed, we will assume that Cd yn is a monotone increasing function of D*. 
It follows from Eqs. (gj) and (J7J) that 

^CL n (x) = G(A n (x), D n (x)) - G(A n -i(x), B n ^(x)). (9) 

According to Eqs. (gj) and (0), the cycle length also satisfies an algebraic condition 

CL n (x) = F(A n _ 1 (x), D n _i(x)) + D n (x), (10) 

and thus Eqs. Q and (fTUjl imply that 

^ [F(A n _x(x), D ft _!(z)) + D n (x)\ = G(A n (x), D n (x)) - G(A n ^(x), D n _x(x)). (11) 

The dynamic pacing protocol gives the following boundary condition at x = 0: 

D n (0) = B- F(A n _i(0), D n _!(0)). (12) 

The sequence of equations Eq. (fTT|) can be solved iteratively to construct Fig. E| If the 
vectors of functions A n (x) and D n _i(x) are known, we can solve (|TT|) to determine D n (x). 
Note that A n+ i(x) can then be computed by applying Eq. ©. 
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FIG. 6: Deflection of the (n + l) st wavefront and waveback due to a premature (B\ < B) S2 
stimulus. The shortened diastolic interval at the stimulus site slows the propagation speed. Solid 
curves represent wavefronts and dashed curves represent wavebacks. 

B. S1-S2 pacing protocol 

In the S1-S2 protocol, tissue is paced at a basic cycle length B until steady-state is 
reached. Then, an 5*2 stimulus is introduced at an interval B\ = B ± 5 following the last SI 
stimulus and the response to the S2 stimulus is measured. In what follows, we assume that 
the S2 stimulus is applied prematurely (B\ = B — 5) following a train of n SI stimuli. The 
S2 stimulus causes a deflection in the {n + l) st wavefront and waveback as shown in Fig. El 

These assumptions imply that 



and the boundary condition 
Equation (fTTj) reduces to 



A n (x) — A , 

D n (x) = (D n (x),D*,...,D*), 
D n (0) =B X - A*. 



(13) 



(14) 



dD n {x) 

dx 



G{A*,B n {x))-G*, 



(15) 



where G* = G(A*, D*) = c d ^ n . Linearizing Eq. (JTiJ) about the point (A*, D*), we have 



where 



dP n (x) 
dx 



A 



-\{D n {x)-D*), 



dG 
11 



(16) 



(17) 



(A*,D*) 



Since we assumed that the dynamic velocity RC is monotone increasing, it follows that 
A > 0. The solution of the Eq. (fTBj) with the boundary condition (JTljl is 

Djx) — D* — 5e~ Xx . (18) 

Let cj^ 2 nt and c b ^ k denote the wavefront and waveback velocities of the action potential 
generated by the S2 stimulus. In order to compute c^ 2 nt , observe that (see Fig. [HJ) 

(j) n+1 (x) = (3 n (x) + D n (x). (19) 

We know that (f> n (x) and /3 n (x) are parallel since they represent the wavefront and waveback 
associated with the final SI stimulus. Therefore, d(3 n /dx = dcf) n /dx = G*, and differentiating 
Eq. (|T9"j) with respect to x gives 

4o„t = G * + SXe -x x - ( 2 °) 
Similarly, to determine c h ™k, we use the expression (see Fig. EJ) 

Pn +1 (x) = 0„+i(x) + A n+l (x) = (j) n+ i(x) + F(A n (x), D n (x)). (21) 

According to Eq. (|13p. the only non-constant argument of the function F is D n (x). There- 
fore, differentiating Eq. (|21|) with respect to x gives 

d/Ul _ #n +1 + |F ^ = ^ + 5Ae _ Ax / x + dF\ ^ (22) 



(23) 



cfe dx dD n dx \ dD n 

which implies that 



c 



Sl2 

back / 

G* + 5\e-* 1 1 + |£ 



The partial derivative dF/dD n in Eq. (J2HJ) is evaluated at (A*,D* - 5e~ Xx , D*, . . . D*). 
Equations (j2*Uj) and (}2*3"|) are analytical expressions for wavefront and waveback velocity for 
the S1-S2 pacing protocol. 

Both formulas (|2*0|) and (}2*3*j) require that we know formulas for the dynamic velocity RC 
(since G* = c~j" ) and the function F. The only difference between the two formulas is the 
presence of the multiplier ^1 + -§§^\ ■ In simple mapping models for which Eq. (J2J) applies, 
the partial derivative 8F/dD n in Eq. (|23|) is replaced by a total derivative F'(D n ). Note 
that the wavefront and waveback velocities approach because A > 0. If the 

S2 stimulus is premature (£>i < B), formulas (|2*0|) and (|2*3*|l show that c b ^ < c^ 2 nt and 
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the pulse broadens as it propagates. Likewise, if the S2 stimulus is late (B\ > B), then 
c back > c front an d the P mse contracts as it propagates. If B\ = B, the formulas reduce to 
c back = c front = c dyn as one would expect. 

Equations (|2(Jj) and (J23)) reinforce our main point: APD restitution, not memory, is 
responsible for velocity rate-dependence. Regardless of how much memory is included in 
the mapping model, Eqs. (J2TH) and (H3J) depend upon D n (x) and no other preceding states. 
The partial derivative dFldD n in Eq. (}2"3*|) represents the slope S12 of the S1-S2 APD 
RC as demonstrated in j9(. As DI decreases, S12 typically increases, thereby increasing 
the discrepancy between the wavefront and waveback velocities. It follows that, in the 
absence of wavefront velocity rate-dependence, APD restitution leads to waveback velocity 
rate-dependence . 

C. An example: Rate- dependent velocity and the two-current model 

In this Subsection, we explain how to apply Eqs. (f2"U|) and (|2*3*|) . using the two-current 
model as an example. As mentioned above, Eqs. (}2T)j) and (|23|) require that we provide 
formulas for the function F and the dynamic velocity RC. Leading-order expressions for F 
and 

Cdyn Can be derived analytically for the two-current model (see Appendix). 

The dynamic velocity RC is provided by Eq. (jA.9|) . Combining Eqs. (J2(J|) and (jA.9|) . we 
generate all S1-S2 wavefront velocity RCs. Likewise, combining Eqs. (|2^|). (|A.3J) . and (|A.9|) 
allows us to construct all S1-S2 waveback velocity RCs. 

All of the analytically derived RCs are shown in Fig. [7| Figure 1?^ shows all wavefront 
velocity RCs. The dynamic and S1-S2 wavefront velocity RCs are indistinguishable. The 
waveback velocity RCs are shown in Figure [7)d. Note the presence of rate-dependence, as 
evidenced by the splitting of the dynamic (solid) and S1-S2 (dashed) RCs. We remark that 
Fig. [7| shows excellent quantitative agreement with the results of numerical simulations 
shown in Fig. 01 

IV. CONCLUSIONS 

We have demonstrated that rate-dependent waveback velocity restitution can exist even 
in the absence of memory. Our numerical simulations show that both the two and three- 
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FIG. 7: Wavefront and waveback velocity RCs generated using Eqs. {2D|! , (ESI), <TO> - and (fOJ) . 

(A = 40 ms, 8 = ±20 ms). Velocities were measured at x = 2.5 cm. (a) Wavefront velocity RCs: 
the dynamic and S1-S2 curves appear to coincide, (b) Waveback velocity RCs: the dynamic curve 
is solid and local S1-S2 curves are dashed. 

current models exhibit rate-dependent waveback velocity, whereas neither model exhibits 
rate-dependent wavefront velocity. We offer a mathematical explanation for the differences 
between wavefront and waveback dynamics. Specifically, comparison of Eqs. (J2(J|) and (|23j) 
shows that as the slope Si 2 of the S1-S2 APD RC increases, the difference between the 
wavefront and waveback velocities is magnified. Therefore, if the S1-S2 wavefront velocity 
RCs coincide with the dynamic velocity RC, then the S1-S2 waveback velocity RCs cannot. 
Moreover, Eqs. (j2"Uj) and (J23J) predict that the splitting between the dynamic and S1-S2 
waveback velocity RCs should be most pronounced at small DI values where S12 is largest. 
These analytical predictions are consistent with the results of our numerical experiments. 
Finally, the validity of the computations in Sec. Ill is strongly supported by the quantitative 
agreement between numerical and analytical investigations of the two-current model (Figs. 
HJandEJ). 
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APPENDIX 

A detailed analysis of the two-current model equations appears in [[18]. With the two- 
current model, Eq. (JTJ) reads 

dv d 2 v h 2 .„ . v . . 

dh ie «<«- (A . 2) 



dt 



Teh 



— v > v ciit , 



where v is transmembrane voltage (scaled to range between and 1) and h is a gate variable. 
The parameters r; n , r c i ose , r out , and r open are time constants associated with different phases 
of the action potential. The gate opens or closes according to whether v exceeds the threshold 
voltage v cr it- Typical choices for the time constants and critical voltage are: r- m = 0.1 ms, 
r out = 2.4 ms, r open = 130 ms, r close = 150 ms, and v crit = 0.13. 

A leading-order estimate of the APD RC is derived in [18j. If the time constants satisfy 
an asymptotic condition r in <^ r out <^ r open , r close , then 

A n+l = F(D n ) = r close In (^^) (A.3) 

\ ^min / 

to leading order, where 

Dn 

h s (D n ) = 1 - (1 - h min ) e -p-, (A.4) 

and hum = 4r in /r out . 



To derive a leading-order estimate of Cd yn , we follow Murray 27]. Assume the fiber 
is paced at a constant basic cycle length until steady-state is reached so that all pulses 
propagate with speed Cd yn = Cd yn (-D*)- We seek traveling wavetrain solutions to Eq. (|A.1J) . 
In the neighborhood of a wavefront, introduce the coordinate 

e = — (*+— V (a.5) 



7"in \ Cdyn 

15 



where the speed Cd yn is to be determined. Assume that v(x,t) = V(£) and h(x,t) = H(£). 
Since r- m is small relative to the time constants in Eq. (jA.2|) . we may safely approximate 
the value of h by a constant in the narrow wavefront region: h « h* = h s (D*). Inserting 
v(x,t) = V(£) into Eq. (jA.l|) . we obtain an ordinary differential equation 

_ V ' + h* V (V. -V){V- V+) = 0, (A.6) 

C dyn T in 

where primes denote differentiation with respect to £ and 




(A.7) 



We remark that V_ is an unstable equilibrium of Eq. (jA.6j) corresponding to the threshold 
for excitation, and V+ is an unstable equilibrium associated with the excited state. We seek 
solutions to Eq. (|A.6|) such that V(£) — > as £ — > — oo and V(£) — > V + as ^ — > oo. It is 
possible to find a solution of a simpler differential equation 

V' = -aV (V - V+) (A.8) 

that also satisfies Eq. ()A.6jl for unique values of the constant a and the speed Cd yn - Substi- 
tuting (|A.8jl into Eq. ()A.6jl . one finds that 



Cdyn 




(A.9) 
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